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(X) ' Abstract 

o ■ 

We perform a time-resolved spectral analysis of bright, long Gamma-ray burst 
; GRB 061007 using Suzaku/WAM and Swift/BAT. Thanks to the large effective 

^ \ area of the WAM, we can investigate the time evolution of the spectral peak en- 

ergy, -Epeak ^^^d the luminosity L\^^ with 1-sec time resolution, and we find that the 
time-resolved pulses also satisfy the E^^ak — Lisa relation, which was found for the 
time-averaged spectra of other bursts, suggesting the same physical conditions in 
each pulse. Furthermore, the initial rising phase of each pulse could be an outlier of 
this relation with higher -Epe^k value by about factor 2. This difference could suggest 
that the fireball radius expands by a factor of 2 — 4 and/or bulk Lorentz factor of 
the fireball is decelerated by a factor of ~ 4 during the initial phase, providing a new 
probe of the fireball dynamics in real time. 
Key words: Gamma-ray bursts 
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1. Introduction 



Many characteristics of Gamma-ray Bursts (GRBs) have been revealed by the previous 
observations, and now it is widely believed that GRBs are the one of the most powerful explosion 
in the Universe, originating at cosmological distances. The association with the energetic 
supernovae was found for some long duration GRBs. This provides a strong evidence that 
the progenitor of the long GRBs is the core collapse of a massive star. The observed prompt 
gamma-ray spectra are often fitted well by a smoothly connected broken power-law called as 
the Band function (Band et al. 1993). However, we are still not able to uniquely associate the 
measured properties to the radiation mechanism of the prompt gamma-ray emission. 

Although, the observed behavior of GRBs, such as the intensity and the peak energy of 
the uFy spectrum (-Epeak) are dramatically different from burst to burst, some clear correlations 
between these parameters were reported, such as the correlation between isotropic equivalent 
radiation energy Eif^o and the peak energy in the rest frame of the burst (-Epcak.src) (Amati et 
al. 2002; Amati 2006), or the isotropic equivalent peak luminosity Lpeafc.iso and the -Epcak.src 
(Yonetoku et al. 2004), and the total radiation energy corrected for the jet opening angle and 
the -Epoak.src (Ghirlauda et al. 2004b). These relations are often used to estimate the inferred 
redshift value, and Yonetoku et al. (2003) estimated redshift for large BATSE sample and 
found that some of them have large redshift with z > 10. Such high-z GRBs are important 
to constrain the cosmological parameters. Furthermore, the same kind of correlation have 
been reported not only for each burst but also for individual pulses of a burst internally. 
Several authors studied this hardness-intensity correlation (HIC) for time-resolved spectral 
analysis. Golenetski et al. 1983 found a power-law correlation between observed energy flux 
and the spectral characteristic energy Eq, and the power-law correlation index was found to be 
a typical value of 1.5-1.7. Kargatis et al. (1994) confirmed this correlation but the power-law 
index showed a larger spread. They obtained the power-law index of 2.2±1.0. Strohmayer et 
al. (1998) also found that this power-law correlation is valid using Ginga data. Borgonovo & 
Ryde (2001) also studied this correlation by strong BATSE sample. They found the specific 
feature that there is the sharp transition from one power-law track to the second track with 
the same power-law index in the decay phase of the pulse. They explained this track-jump 
behavior as the transition to a new pulse. The fact that the HIC holds for both each pulse and 
within the pulse implies that the same physical condition is realized in each pulse. Therefore, 
to investigate this relation it is very important to consider the radiation mechanism of the GRB 
prompt emission. Recently, Liang et al. (2004) analyzed the Liso and -Epcak.src of the 2408 time- 
resolved spectra for 91 BATSE GRBs, assuming that the burst rate as a function of redshift is 
proportional to the star formation rate. Yoshida et al. (2007) measured the Liso and -Epeak.src 
value for time- resolved spectra of HETE-2 bursts. Both of them found that L-^so oc -E'pcak.sro 
and Yoshida et al. (2007) found that this relation is tighter for each spike-averaged spectra 
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rather than that of simply 5.2 s-divided spectra. Therefore, this relation could be different 
for finer time-resolved spectra. However, most of previous analysis was done using data for 
each pulse or for the decay phase of such pulses, while the detailed properties of time-resolved 
spectra, including initial rising phase for each pulse have not been considered and thus are still 
ambiguous. In this paper, we report on the time-resolved spectroscopy of the bright long GRB 
061007 observed by the Suzaku/WAM and Swift/BAT. This burst was also observed by many 
ground-based telescopes and the redshift was measured to be 1.261 (Osip et al. 2006; Jakobsson 
et al. 2006), and Schady et al. (2007) suggested that this burst has an early jet break time 
which is within 80 s of the prompt emission, and a highly collimated outflow is likely case for 
this burst. Using this bright burst with a measured redshift, we can investigate the relation 
between Ljso and -Epcak,src within each pulse, including both the pulse-rise and the decay phase. 
This is possible since Suzaku/WAM has the largest effective area from 300 keV to 5 MeV than 
any previous missions. Furthermore, since Swift/BAT covers an energy range from 15 keV to 
150 keV, the joint analysis using Suzaku/WAM and Swift/BAT gives us further constraints 
on the spectral parameters. Therefore, this joint, broad-band Suzaku/WAM and Swift/BAT 
analysis allows the excellent time-resolved spectral analysis with both fine time resolution and 
high sensitivity. We also investigate the time-averaged properties of this burst such as E^go - 
-E'peak,src (Amati) relation and L^so - -Epeak.src (Yonetoku) relation in order to compare with other 
bursts. Thus, we denote our "time-resolved" peak energy and luminosity relation as -Epeak.src ^ 
L\gg relation in this paper to distinguish from those time-averaged relationships. 

2. Observations 

2.1. Suzaku/WAM 

The Suzaku/WAM (WAM) is the active shield of the Hard X-ray detector (HXD-II) 
(Takahashi et al. 2007, Kokubun et al. 2007) onboard Suzaku (Mitsuda et al. 2007). It consists 
of large and thick BGO crystals and it is also designed to monitor the all-sky flux over the 50 keV 
to 5 MeV band. It has the largest effective area from 300 keV to 5 MeV of all GRB missions 
and thus enables us to perform a wide-band spectroscopy of GRBs with a high sensitivity 
(Yamaoka et al. 2005, Yamaoka et al. 2006). There is still uncertainty in the detector response 
matrices of the WAM and the observed flux fluctuates by 10-20% on average and about 40% at 
the maximum depending on the incident azimuthal (0) and zenith (6) angle of GRBs (Ohno et 
al. 2005). We found that this uncertainty (as compared to the pre-fiight calibration) is mainly 
caused by the absorption by the structures attached inside the satellite panel. We confirmed 
that this trend is the same even in the in-orbit environment by the cross-calibration between 
Swift/BAT and Konus-Wind, using the data of simultaneously detected GRBs (Sakamoto et 
al. 2008 in preparation). 

The WAM outputs two data types, the transient (TRN) data and the gamma-ray burst 



3 



(GRB) data. The TRN data are always accumulated with 1 sec time resolution and 55 energy 
channels. This can be used to monitor the bright soft gamma-ray sources with the Earth 
occultation method, as was done by the CGRO/BATSE (Ling et al. 2000). On the other hand, 
the GRB data are recorded for 64 sec only when the GRB trigger is activated, and the data 
cover 8 sec before and 56 sec after the trigger time. The GRB data have four energy channels 
with 1/64 sec time resolution, in addition to the spectral data in 55 pulse height channels with 
0.5 sec time resolution. 

GRB061007 triggered the WAM at 10:08:05 UT, October 07, 2006 (To(WAM)). Three 
large peaks as well as a following smaller peak were clearly seen in the WAM light curve as 
shown in figure 1. Gamma-ray photons are strongly detected by the WAM-3 detector. The Tgo 
duration of this burst was about 59 seconds in the 50-5000 keV band (Yamaoka et al. 2006b). 
Since this burst was so long that the emission lasted more than 56 s after the WAM trigger, the 
GRB data system could not record the whole emission profile. This is why the light curve of the 
WAM GRB data look truncated above 56 s (see Figure 2). Based on the position information 
of this burst provided by the Swift as mentioned below, the incident direction against the WAM 
detector was found to be (^, 0) = (92°, 198°). This incident direction is consistent with the 
strong detection by the WAM-3 detector, and the detector response of this direction is known to 
be reliable by the pre-flight and in-orbit calibrations. The uncertainties of the response matrix 
should be within 20-30% for this burst. 

2.2. Swift/ BAT 

The Swift/BAT (BAT) was also triggered by GRB 061007, 3 s after the WAM trigger, at 
10:08:08 UT, October 07, 2006. The refined BAT position for this burst was RA = 03''05'^llf8 
Dec = -50°29'47."7 (J2000) with an uncertainty of 0.9' at 90% containment (Markwardt et 
al. 2006). The BAT light curve showed a multi-peaked structure with a duration of about 100 
seconds. The Swift/XRT and Swift/UVOT observed this burst 80 s and 400 s after the BAT 
trigger, respectively, and they found a bright fading source. 

3. Data Analysis and Results 

Although, the main purpose of this study is to investigate the behavior of the time- 
resolved spectra of GRB 061007, before we perform the time-resolved spectral analysis, the 
time-averaged spectral parameters have been also investigated. We utilize the transient (TRN) 
data for the WAM because GRB data of the WAM do not contain the whole emission of this 
burst due to its long duration. Therefore, we can analyze only 1-sec time resolution data for 
this analysis. We perform the spectral analysis by applying three spectral models. The first 
one is a simple power law (PL) model: 




(1) 
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where A is the normahzation constant at 100 keV in photons cm~^ s~^ keV and a is the 
power law photon index. 

The second model is a power law with an exponential cutoff (CPL): 

where -Epcak is the peak energy in the vFy spectrum and it represents the energy at which most 
of the power is emitted. 

The third model is a smoothly connected broken power law model known as the Band 
model (Band et al. 1993). 



N{E) = A X 



VlOOkeVy/ ^Peak , 



(2 + a) 

/ K 

A X 



VlOOkeV; 



^ (a -/9)-gpcak 

100keV(2 + a) 



exp(/5 — a), 



for E > ^^7^^^r' (3) 

where a is the power law photon index in the lower energy band, and /5 is that in higher energy 
band. In addition, in these basic models, we add the constant factor when we perform the joint 
spectral analysis between the WAM and the BAT to take into account for the uncertainties of 
the detector response matrix of the WAM. Throughout this analysis, we fit the spectrum from 
120 keV to 5000 keV for the WAM and from 13 keV to 150 keV for the BAT, respectively. 

In order to extract the dead-time corrected light curve and the spectrum, we use the 
standard FTOOLS; hxdmkwamspec and hxdmkwamlc for transient data, and hxdmkbstlc and 
hxdmkbstspec for GRB data, which are included in the HEAsoft software package version 
6.0.6 (http://heasarc.gsfc.nasa.gov/lheasoft/). We also apply a 2% systematic error to the 
WAM spectrum only for the time-averaged spectral analysis due to its high counting rate. We 
also derive the light curve and the spectra from Swift/BAT data for joint spectral analysis. 
When we extract the spectrum or light curve from the BAT event data, we use the standard 
FTOOLS; batbinevt, after we re-processed the energy calibration, making the detector quality 
map and mask weighting, using the latest version of the software and calibration data base. We 
also apply systematic error for the BAT spectra, using batphasyserr script. The systematic 
error vector is retrieved from the calibration data base. 

We use the XSPEC version 11.3.2 for spectral analysis (Arnaud 1996). All quoted error 
are 90% confidence level. 
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3.1. Time- averaged Analysis 

We can see that the emission profile is well separated into first weak FRED like structure 
and second bright, multi-peaked episode. Therefore, firstly, we divide this burst into two 
time regions. The first time region named as interval A is selected from To(WAM)-2.0 s to 
To(WAM)-|-18.0 s, where To(WAM) means the WAM triggered time, and the second time 
region of the interval B is from To(WAM)+27.0 sec to To(WAM)+87.0 sec. We also analyze 
the integrated time region from To(WAM)-2.0 sec to To(WAM)+85.0 sec. We also extract the 
BAT spectra from the same time interval as that of the WAM in order to perform the joint 
spectral analysis. The selected time region is shown in the light curve of the WAM and BAT 
in figure 2. The best-fit parameters obtained by the joint analysis for each time region are 
shown in Table 1. Figure 2 shows the spectra of the joint analysis with the WAM and the 
BAT for whole time region (A+B). We find that the spectrum of this burst can be described 
by the typical Band model for each time interval. We obtain the photon index of low-energy 
part, a of -0.91(±0.07), -0.73(±0.03), -0.79(±0.03), photon index of high-energy part, (3 of 
-3.001?:^^, -3.22l{]i^,-3.2llJ]JJ], and the peak energy, Epeak of 354l| keV, 468 (±10) keV, 460 
(±10) keV, for interval A, B, and A±B, respectively. These spectral parameters are consistent 
with that of reported by Konus-Wind (Golenetskii et al. 2006). The 100-1000 keV fluence and 
the peak flux in 1-s time scale is measured to be 1.6 (±0.15) xlO~^ erg cm~^ and 14.5 (±0.9) 
photons s~^ cm~^, respectively. The constant factor of the WAM against the BAT is about 
1.10 for any time regions. These values are consistent with the current uncertainties of the 
response matrix of the WAM. Furthermore, we find that the spectral parameters can be tightly 
constrained by the joint analysis between the WAM and the BAT. 

We also investigate other spectral properties of this burst using the spectrum of the 
whole time interval (A+B). First, we measure the hardness ratio between 100-300 keV and 
50-150 keV, using the fluence ratio based on the spectral fltting result and we obtain the value 
of 3.52 (±0.13). We plot this value on the hardness- duration plane in Figure 3 (left panel). We 
also plot the BATSE results together for comparison, and we flnd that this burst locate at the 
standard long/soft regime like BATSE results. Then, we investigate the total emitting energy 
of this burst. We calculate the isotropic equivalent total energy Ei^o from 1.0 to 10000 keV 
and the peak energy at the rest frame -Epcak.src, using cosmological parameters (Hq, I^a, ^m) = 
(65, 0.7, 0.3), and we obtain the Ei^o of 1.0(±0.3) x 10^^ g^g, and the ^pcak,src of 1040li keV, 
respectively. We plot this result in the -Epeak,src--E'iso plane and we show it in Figure 3 right. For 
long GRBs, it is well known that there is a strong correlation between Eiso and -Epcak.src (Amati 
relation; Amati et al. 2002, Amati 2006). Thus, we also plot the previous results reported by 
Amati 2006 in the same flgure, and we flnd that the Eiso and -Epcak.src value of GRB 061007 is 
quite consistent with Amati relation. Figure 4 also shows a comparison of the -Eiso and 1-sec 
peak luminosity (Yonetoku) relation (Yonetoku et al. 2004), and we flnd that GRB 061007 also 
satisfles the Yonetoku relation. From those time-averaged spectral analysis, we conclude that 
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there is no large systematic problem in our analysis and that GRB 061007 is a typical long 
GRB with many common characteristics such as the duration, the spectral hardness, and the 
total emitting energy. 

3.2. Time-resolved Analysis 

After we confirmed that there was no problem in the time-averaged spectral analysis 
between the WAM and the BAT, we investigate the time-resolved spectral properties of this 
burst. In this analysis, we divide the WAM data in 1.0 s time resolution. We also divided the 
BAT data for the same time region. We extract total 58 time-resolved spectra and perform 
the spectral analysis with the same manner as that we used in the time-averaged analysis. The 
time-resolved spectra are also described by the GPL or Band model. We can constrain the 
high energy photon index f] from only a few spectra which belong to the bright time region. 
Therefore, we applied the fixed value of f] of —3.2 which is obtained at the time-averaged 
spectral analysis. The best-fit parameters for 58-time-resolved spectra are shown in Table 2 
and 3. From these time-resolved analysis, we find the evolution of the spectral shape during 
the burst. Figure 5 shows examples of the uFi, spectra with the best-fit Band model for various 
intensity region, and the change of the spectral shape around the -Epeak is clearly seen. In figure 
6, we present the best-fit parameters obtained by the Band model fit as a function of time, 
and both the low energy photon index a and the Epcak show a hard-to-soft tracking behavior 
correlated with the burst intensity. The low energy photon index a changes between —1.5 to 
—0.5 with an average value of —0.84, which is consistent with the time integrated parameter 
of — 0.79(±0.03). The Ep^ak also evolves within the burst, and it moves from 120 keV to 633 
keV with an average value of 375 keV. Figure 7 shows the relation between low energy photon 
index, a and the -Epeak, and we find that only one time-resolved spectrum have the low energy 
photon index of flatter than —2/3, which is not allowed by the standard synchrotron shock 
scenario. 

4. Discussion 

4-1- -E'peak ^-^iso i^^litiou foT Ume-resolved spectra 

From time-resolved spectral analysis, we find that the time-resolved -Epcak! -^'pcak changes 
are correlated with the burst intensity. In order to confirm this, we perform a correlation 
analysis between -Epeak and time resolved luminosity; Liso for time-resolved spectra. The L\g^ 
value is measured as the 1-sec time-averaged value for each of the time-resolved spectra in 30 
keV to 10000 keV energy range. Figure 8 shows the relation between -Epeak and L\^^. We can 
see a clear positive correlation between these two parameters. To quantify this, we calculate 
the Spearman rank-order coefficients (r^), and we obtain the correlation coefficient of 0.78. 
This coefficient corresponds to the significance probability that there is no correlation of 5.4 x 
10"^"^. Therefore, we confirm that the correlation surely exists at > 3cr confidence level. We 

7 



then fit this correlation by the simple power-law model like the Hardness-Intensity correlation, 
the Amati relation, and the Yonetoku relation. The fit gives a best fit parameter as, 

/ J^t \ / 7-t \0.43(±0.03) 

= 492(±24) X -^ps^ r . (4) 

\1 keVy ^ ^ V1052 erg s-^ J ^ ^ 

The power-law index slightly flatter than the value obtained by Golenetskii et al. 1983 
but roughly agrees with the previous analysis by Yonetoku 2004; -Epeak oc L?^. Furthermore, 
our data set can constrain both the power-law index and the normalization factor as well as the 
previous analysis using only one burst data set. The value for this fit is 83/56 and we still 
see a certain deviation in this relation around Lj^^/lO^^ erg s~^ = 2 — 5 (see Figure 8), and the 
data point which have large dispersion from this relation tend to have a higher E^^^^ value by 
a factor of 2 compared with many other data. We calculate the E^^^^^f {L\^^)^-'^^ for each time 
region in order to separete these data point from other data which follow the equation(4). If the 
correlation follows -Epcak -^Lo'^^ f*^^ time region, this parameter should retain constant. 
Figure 9 shows the E^^^]^/ {Ll^^)^'"^^ parameter as a function of time. We also show the best-fit 
constant value with 3 cr confidence region. The 8 data points (T=4, 5, 6, 7, 31, 32,48, 50) exceed 
3 a limit of the constant value, and thus we consider these data points as the outliers on 
the -Epeak ~ -^iso pl^ne. We confirm that these outliers with large dispersion from the main 
population are not caused by the systematic uncertainties in the analysis procedures such 
as the uncertainties of detector response matrix, systematic effect to the obtained spectral 
parameters, and the effect of the intense spectral evolution during the initial rising phase. The 
details of this verification processes are shown in the Appendix. 

We then divide the derived parameter pairs into the main population, where the 
-^peak/(-^iso)°'^^ correlation is acceptable, and the outlier population, where that relation is 
broken at 3 a significance. We calculate the correlation coefficient and perform the power-law 
fit for each population separately, and obtain the coefficients of 0.90 (P^-s = 6.9 xlO~^^) for the 
main population, and 0.62 for the outliers. The correlation coefficient of the main population 
is marginally improved when we exclude the outliers. The outliers also have a certain degree 
of the correlation. The power-law best-fit parameters for those two populations are found as 
follows, 

/^t \ / \0.46(±0.03) 

= 456(±25) X for main population 

yl keV / ylO^^ erg s~^/ 

/ x0.25(±0.14) 

= 909(±147) X [j^E^^^,) for outliers. (5) 

For the main population, the best-fit parameter of the power-law fit does not change 
significantly from that we derived from the entire data. We find that the power-law model is 
still valid for this relation. Furthermore, this sparation also led to a very interesting result. 
From figure 9, we can see that all of these outliers correspond to the rising phase of both 
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initial weak peak and the second and third intense pulses. This indicates that if we exclude 
the initial rising phase of each pulse, we obtain much tighter -E'p<,ak,src^-^iso correlation for all 
GRB. Therefore, it is probably more appropriate to exclude initial rising phase of each pulse 
when we use this relation as the redshift indicator like time-averaged -Epeak,src ^ -^iso (Yonetoku) 
relation. Moreover, it is likely that these outliers, in other words initial rising phase of each pulse 
have a different correlation with higher normalization and possibly smaller power-law index. 
However, we cannot give a strong conclusion for this because this power-law index varies, 
strongly depending on the criteria between the main population and outliers. Furthermore, we 
cannot constrain the power-law index of outliers due to the small number of data. Therefore, the 
robust conclusion is that the initial rising phase of burst have a different -Epc^k ~ -^Lo correlation 
with higher -Ep^ak value from that of decay phase of burst. Note that this argumentation do not 
form a circular logic because the main point of this study is not to improve the fit of -Epeak ~ -^Lo 
correlation by removing the outliers but to just find that these outliers correspond to the initial 
phase of each pulse. This kind of trend that the initial phase of long GRBs has harder spectrum 
than that of decay phase is already reported by Ghirlanda et al. 2004. However, the WAM data 
revealed that this harder-rising phase of long GRBs also follow the another power-law relation 
with that of softer-decay phase. This difference between these different phases affects mainly 
the normalization of the relation. This fact can be useful to constrain the emission mechanism 
or dynamics of the emission site of the prompt emission of GRBs. 

4-2. Implication for the emission radius and the bulk Lorentz factor 

Let us consider theoretical implications of our results that the spectral peak energy -E'peak 
evolves from hard to soft in the rising phase of each pulse. Although the emission mechanism 
of prompt GRBs is still enigmatic, the leading models are (1) the synchrotron shock model 
(Zhang & Meszaros 2002) and (2) the photosphere model (Meszaros & Rees 2000; Thompson 
et al. 2007; loka et al. 2007). 

In the synchrotron shock model, the peak energy is identified with the typical syn- 
chrotron energy -E'peak = Th'j'^eB / rrieC of electrons that are shocked in the relativistic outflow 
with a Lorentz factor F. Assuming that a fraction eg of total energy goes into the electron 
acceleration and a fraction goes into the magnetic amplification, we have the typical Lorentz 
factor of electrons as 7m ~ ^emp/nie and the magnetic luminosity as esLtot = 47rr^cF^i?^/87r 
where r is the shock radius and L^ot is the total luminosity. Since almost all electron energy 
is radiated, the photon luminosity is given by = eeLtot- Combining these relations, we have 
the peak energy as a function of r, F and as 

^peak ~ 3 e^J^el'^L^r^i MeV, (6) 

where Qx = <5/10^' in the cgs unit. In this model, the low energy photon index should not 
exceed —2/3, and even —3/2 if electrons cool fast as expected in GRBs (Ghisellini et al. 2000). 
However, many previous observations found flatter low energy index than —2/3 (e.g., Preece 
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et al. 1998; Ghirlanda et al. 2002). Since our results for GRB 061007 exhibit flatter low energy 
photon index than —2/3 for the spectrum of one time interval, the photosphere model should 
also be considered. 

In the photosphere model, the peak energy is identified with the thermal peak -Epeak ~ ^T' 
of the photosphere under which the outflow energy is internally dissipated and thermalized. 
Most opacity of the photosphere is provided by pairs or electrons associated with baryons 
(e.g., Murase & loka 2008). This model has several advantages that the peak energy is stabilized 
even for high radiative efficiency and that the low energy photon index can be as hard as the 
thermal one, ~ 1 (Thompson et al. 2007; loka et al. 2007). With the Stefan-Boltzmann formula, 
we have the comoving energy density of photons as aT"^ = L^/Anr'^cT'^, so that 

£;peak ~ 1 Ty\-,^'^L% MeV. (7) 

From equations (6) and (7), the high peak energy above the Yonetoku relation Ep^ak oc 
Llj"^ suggests that the shock radius r is small and/or the Lorentz factor F is large. This means 
that the fireball is expanding and/or decelerating in the initial rising of pulses, and we could be 
directly observing the fireball dynamics through the evolution of the Yonetoku relation. Since 
the shift of the peak energy is about a factor of two, the radius increases by a factor of 2-4 
and/or the Lorentz factor decreases by a factor of ~ 4. Because the time dependence could be 
different between the synchrotron shock model E^cak oc oc and the photosphere model 
Epeak F^/^r~^/^ oc F^/^t~^/^, we could potentially discriminate models with similar but more 
detail spectral observations. 
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Appendix 

In this section, we present the detail of the verification processes to check whether these 
outliers in the Figure 8 are real or not from the point of view of the spectral analysis. At 
first, we compare the spectral shape of these outliers with that of main population as shown 
in figure 11 left. For this comparison, we produce two integrated spectra. One of them is the 
spectrum of outliers which are summed up among all 8 data points which we defined as the 
outliers, and the other is the spectrum of the main population which is the summed up 9 data 
points with similar luminosity as the outliers. These selected data points are shown in the same 
figure. From this comparison, we find that the spectral shape of outliers is surely different from 
that of main population, and thus we confirm that these outliers do not simply arise from the 
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uncertainty of the fitting procedures. Then, we investigate another possibihty that the lack of 
time resolution of the WAM-TRN data causes the behavior that the pulse rising phase tends to 
have higher ^'peak value. We can only resolve with 1-sec time scale by the WAM data even if the 
hard-to-soft spectral evolution exists with finer time scale. The -Epeak value might be shifted if 
we average the spectrum out with 1-sec time scale which have that hard-to-soft evolution. In 
order to confirm this, we performed a Monte Carlo simulation to reproduce the spectra with 
hard-to-soft evolution. In this simulation, we assume that the E^^^y. and the luminosity of each 
spectrum change, following -Epcak oc L?,^, as in the Yonetoku relation, to represent the spectral 
evolution. We then investigate the differences of the -Ep^ak value between each spectrum and 
that of averaged spectrum. As a result, we find that if we average the spectra with hard-to-soft 
evolution, the -Epeak value becomes somewhat higher than than of each simulated spectrum, by 
about a factor of 1.3. However, for the case of GRB 061007, the difference of the -Epcak value 
between the outliers and the main population is much larger than this effect, which is about a 
factor of 2, and we cannot explain the behavior of the outliers only by this spectral evolution 
effect. 

Next, we estimate the systematic uncertainties of these time-resolved spectral param- 
eters. There are three major systematic uncertainties for the WAM data. The first one is 
the uncertainty of the detector response matrix. However we already confirm that this un- 
certainty should be within 20% from the result of time-averaged spectral analysis. Moreover, 
since we compare each time-resolved spectrum of only one burst, the uncertainties of response 
matrix do not influence our time-resolved analysis. The second issue is the dependence on the 
burst intensity. To investigate this uncertainty, we calculate various types of spectra with same 
spectral parameters and different intensities by Monte Carlo simulation. Figure 12 left shows 
the variation of the spectral parameters as a function of the burst intensity obtained by these 
simulated spectra, and we find that the -Epcak goes to a lower value as the burst intensity goes 
down. However, this variation should be within only 20% and this tendency is inconsistent 
with that of outliers of our time-resolved spectra, i.e. the -Epcak value of outliers becomes larger 
than major population. The third uncertainty is the variability of the background. Figure 12, 
right panel, shows a comparison of the background spectra between pre- and post-burst time 
intervals, and we find that the background spectrum only varies by less than 5%. We find that 
the -Epcak value changes only 10% at maximum even when we use the background larger by 5%. 
This amount of variation cannot explain the differences of the -Epeak between the main popu- 
lation and the outliers of about a factor of 2. From these estimation of the systematic effects 
of our time-resolved spectral analysis, we conclude that the differences of the main population 
and outliers cannot be explained by the systematic uncertainties and this different population 
is real for this burst. 
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Fig. 1. Light curve of GRB 061007 obtained by the WAM GRB data in 1/64 s 
time resolution. This hght curve is divided into each energy band of the GRB 

data (THO: 50-110 keV, THl:llG-240 keV, TH2:240-520 keV, and TH3:52G-5G00 keV). 
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Fig. 2. Top: Light curve of GRB 061007 obtained by the WAM monitor data (TRN data; top) 
and the BAT (bottom). The time regions for the time averaged spectral analysis are also shown 
by the vertical dashed lines. Bottom: Joint spectral fitting with a Band model for the WAM (red) 
and the BAT (black) data for the data A plus B. The solid line represent the best fit Band model. 
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Table 1. Results of the time integrated joint fittings for GRB 061007. 



Model 


Det 


a 13 




C(WAM)1" 




Interval A 


PL 


WAM3 

R AT 

BAT+WAM3 


t-uy_o.o5 




^^•'^^-0.14 


54/24 

/ "J ' 
570/82 


CPL 


WAM3 
BAT 

BAT+WAM3 


1 1O+0.34 
^•^''-0.38 

0.92tro? 


3681:^^ 
366li^^ 


1 07+°-^° 

J^-'J'-0.09 


15/23 
48/56 
65/81 


Band 


WAM3 
BAT 

BAT+WAM3 


n 07+0 22 
_Qqi+0.07 _2 05+0.43 


358i« 
3561^3^ 


l-U»_o.o9 


13/22 
48/55 
62/80 


Interval B 


PL 


WAM3 
BAT 

BAT+WAM3 


2.05 

n oQ+o.oa 

1.71 




2.31 


639/24 
21/57 
4052/82 


CPL 


WAM3 
BAT 

BAT+WAM3 


0.80l«;i^ 
0.92±H^3 

n 75+0.03 

'J- ' "^-0.03 


478ti^ 
4781 


1 19+0.05 

^•-'^^-0.04 


36/23 
21/56 
65/81 


Band 


WAM3 
BAT 

BAT+WAM3 


-n 66+"'^* -3 22+"'2^ 

'-'•""-0.16 'J-^^-0.42 

-u.az_Q Q3 

|-| 7/1+0.03 q qi +0.27 
'-'■'^-0.03 '-'■'^^-0.45 


460lfg 
4671 


1 10+0. 05 
-'-•-'-'^-0.04 


22/22 
21/55 
52/80 


Interval A+B 


PL 


WAM3 
BAT 

BAT+WAM3 


2.04 

1.66 




1.87 


456/24 
27/57 
3169/82 


CPL 


WAM3 
BAT 

BAT+WAM3 


0.83l»;ii 
0.911H? 
0.82iH^ 


4671^6 
46811^ 


1 19+0.05 

^•^^-0.05 


29/23 
27/56 
61/81 


Band 


WAM3 
BAT 

BAT+WAM3 


-0.681°;?° -3.17l°;5« 
-0 qs+^ o^* 

'-'■^°-0.03 

r|Oi+0.03 097+0.29 
'-'■°^-0.03 '-'■^'-0.54 


4471?° 
456i;« 


1 10+0. 05 

J^-J^'J-o.os 


18/22 
26/55 
51/80 



f: constant factor of the WAM against the BAT in joint spectral fittings. 
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10-" 10-3 0.01 0.1 1 10 100 1000 10" 
E,,„(10S2 erg) 

Fig. 3. Left panel shows the relation between Tgp duration and hardness ratio (100-300 keV to 50-100 
keV). The dots show BATSE data (Paeiesas et al. 1999). The value of GRB 061007 is plotted by the 
red cirele. Right panel shows the i?poak,src-£'iso relation of the previous observations (Amati 2006). 
The solid line shows the power-law best fit of this correlation; E'pcak.src = 95 x E^^^^. The dashed 
line shows the 2 a confidence region. The result of GRB 061007 is indicated by the red circle. 



10" 




0.1 1 10 100 

Lp=.K(1sec) [10=2 erg s-'] 

Fig. 4. Relation between E'pcak in the source frame and peak luminosity of 1-sec time scale. 
Previous result obtained by Yonetoku et al. (2004) is shown by black circles, and our re- 
sult of GRB 061007 is shown by red circle. The solid line and dashed-lines show the pow- 
er-law best fit of this relation; Lisa — 2.34 x 10^''' x -Epoaksrc ^^'^ ^ a confidence region. 
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unfolded spectrum 



Swift/BAl' Suzaku/WAM 




channel energy (keV) 

Fig. 5. lyFi, spectra of three time-resolved spectra of GRB 061007. These spectra are extracted 
by 1-sec time resolution. Solid lines represent the best fit Band modes. Black squares, red cir- 
cles, and green diamonds correspond to the data of the time region of T=10, 36, 51, respectively. 
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Table 2. Best-fit parameters of Band model for 58 time resolved spectra of GRB 061007. 



Band model parameters Intrinsic parameters 



Time 


a. 


4cak[keV] 


C(WAM)''' 




-Epcak.src [keV] 


LLo erg s-i] 


2 


-o.o4l?:I^ 


J'J^_209 




79/81 


1209+4°f 




3 


-0.86ir4? 


246t??^ 


1 70+137 

i-'y_0.94 


75/81 


c;t:7+768 
00 < _261 


0.23to.n 


4 




447+373 


1 42+0.88 
J^-^^-0.51 


78/81 


10071^41 


0.68«-J^ 


5 


-0.85l°:i 


596+?« 


1 07+0.44 
^••J ' -0.32 


68/81 


1349lf,g 


1.96«:g 


6 


-o.83ts;:f. 


495i^« 


1 =9+0.36 


88/81 


1119113 


2 t;c+0.53 


7 


-o.89l°:l? 


456^^" 


0.991^19 


71/81 


1031^? 


9 O/I+0.58 
^•'J^_0.46 


8 


r, 79+0.25 
' ^-0.23 


2691^8 


0.97ig:g 


84/81 


609tii 


1 9r+0.36 
1^-^0_Q 27 


9 


-l-08t°il 


265i|f 


i-24t«:f3 


80/81 


599t?^i 


0.85t°J} 


10 


-0.49l«:- 


390lf3 


0.86l«:f, 


89/81 


O09+148 


9 c;9+0.61 
^•0^-0.49 


11 


-0.69l°:l^ 


282+4^ 


0.98l°:?^ 


97/81 


638+^f 


1 76+°-^^ 

^- '"-0.30 


12 


-0 82+°-2'' 


238t^^ 




75/81 


c;qo + 163 


1 n5+°-2^ 

i.ua_g 22 


13 


-0 7S+° i^ 


2511^3 


1 01+0.29 

J^-JJ^_0.23 


72/81 


5691^6 


l-DO-0.26 


14 


1 19+0.19 
^•^^-0.18 


274l«^ 


1 46+°-^^ 

J^-'±'J_0.29 


89/81 


619ti^? 


1 97+0.28 
^•^' -0.22 


15 


-0.70l«j^ 


163i^?« 


^^•^^^-0.37 


67/81 


070+249 

'J'^-130 


0.56t°:?« 


16 


-0.94l«:^| 


2241^? 


J^-J^*-0.30 


63/81 


507t?i^ 


o.solgj^ 


17 


1 in + 0.32 
^•^■^-0.27 


197ljf 


^^•^"-0.38 


85/81 


44515^7 


o.eilg;- 


18 


o.iii?:S 


85t^f 


2-5it?:f7 


82/81 


19ll^r 


22+° " 

"•^^-0.06 


19 


-1 fll +2-23 
^^•"J^-0.72 


ii4t^r 


o-9ot^:i 


97/81 


-140 


iq+0.30 



f: constant factor of the WAM against the BAT in joint spectral fittings. 




800 



Epeak(keV) 



Fig. 7. The relation between observed peak energy i^p^ak ^^id low energy photon index ol. The 
limit of a value predicted by the synchrotron emission model (a = —2/3) is indicated by horizontal 
dashed-line. Upper region above this line is not allowed by the standard synchrotron emission model. 
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Table 3. Continued. 



Band model parameters Intrinsic parameters 



Time 


Q 




4cak[keV] 


C(WAM)''' 


xVd.o.f 


-Epeak.src [keV] 


LLo [10'= erg s-!] 


30 


-1. 


1 1+0.72 
■-'--'--0.47 




-, 07+1.43 

r-o/ _Q g3 


65/81 


7f!4+4148 
' "'1-489 


""^"-0.18 


31 


-0. 


77+0.43 
'-0.36 




1 4q+0.80 


92/81 


= 97+237 
>1^' -159 


"•^"-0.14 


32 


-0. 


64+0-20 




J^-^^-0.23 


92/81 


13541??? 


9 q7+o.67 
^•^ ' -0.55 


33 


-0. 


= 9+0.16 




^^•"'-0.16 


75/81 


1310l??7 


r 90+0.87 


34 


-0. 


■O'3-O.ll 


571IS 


J^-^J^-0.14 


81/81 


1292ljf 


6 49+0-85 

"•^^-0.76 


35 


-0. 


■'^'^-0.10 


408t^^ 


1 97+0.15 
^•^ ' -0.13 


104/81 


9231S 




36 


-0. 


70+0.09 
■ ''^-0.09 


4481^3^ 


1 21+°-" 

^•^^-0.12 


71/81 


10131^5 


r or+0.63 
^•0^-0.57 


37 


-0. 


7-1 +0.08 
■ ' ^-0.08 


6561^^ 




88/81 


1485llJ^ 




38 


-0. 


Rq+0.08 
■"^-0.08 


553 


i.i8i«:- 


74/81 


12501^^ 


oq+0.83 


39 


-0. 


7e+0.0S 
■"^-0.08 


579t« 




80/81 


1310l^t 


8.36l°-;«3 


40 


-0. 


qn+0.09 
■^^"-0.08 


469tf2 


-1 07+0.15 

^••^' -0.14 


69/81 


loeil^f 


4.75l°:f7 


41 


-0. 


7Q+0.08 
■ '^-0.08 


574lf3 


1 90+0.12 
r-^o_Q 11 


84/81 


12991^°^ 


7 9r+0.73 
' -^^-0.66 


42 


-0. 


79+0.07 
■ '^-0.07 


643j:*| 


-1 -r 0+0.10 
^•^■^-0.09 


103/81 


1456l^«* 


9 99+0-91 


43 


-0. 


■"^-0.07 


507^^3^ 


i.28i«:l? 


103/81 


ii47l!^i; 


7 34+0.67 
'"J^-o.ei 


44 


-0. 


74+0.07 

■ '^-0.07 


623lg 


^•^^-0.10 


81/81 


1410l^« 


10 io+°-*'' 


45 


-0. 


70+0.07 
■ '"-0.07 


5871^3^ 


1 25+'^-" 
J^-^J-o.io 


66/81 


1329lg? 


9 70+0-82 
"-'-0.76 


46 


-1. 


■'^"-0.10 


333t^^ 


-1 9r+0.19 


78/81 




i.yu_Q 25 


47 


-1. 


-r r+0.22 
■J^'J-0.18 




1 01 +0.40 
r-Ui_(j 28 


70/81 


596l|g 


0-63lo.?? 


48 


-1. 


07+0.14 
■"'-0.14 


352t|f 


r-^u_Q 23 


72/81 


7971?^? 


1 07+°-2'^ 

l-"'-0.20 


49 


-0. 


04+0.14 
■O4_o.l3 


2941™ 


-1 90+0.24 
-'-•^°-0.20 


86/81 


66211?! 


1 -^4+0-27 
l"-'*-0.20 


50 


-0. 


QO+O.IO 
■'^'^-O.IO 


412l^« 


J^-^^-0.16 


81/81 


qo-|+178 


2.06l°J^ 


51 


-0. 


.0+0.06 
■48-0.06 


563127 


1 29+°-'^° 


62/81 


12741^? 


10.941°:?^ 


52 


-0. 


74+0.07 
■ '^-0.07 


4101^2 


1 09+0-09 


95/81 


9271?^ 


"•ri_o.49 


53 


-0. 


OC-+0.07 
■°'^-0.07 


399I34 


r-^y-o.ii 


52/81 


903l?8 


4 QQ+0-41 
^•'J'J-0.37 


54 


-0. 


70+0.07 
■ '''-0.07 


009+32 


1 97+0.11 
^•^ ' -0.10 


60/81 


8651^^ 


*-"O-0.37 


55 


-0. 


00+0. 06 
■00-0.06 




-1 -10+0.10 

^•^■^-0.09 


83/81 


7981?? 


4 59+0--11 

^•"^-0.37 


56 


-1. 


Qg+0.09 
■"'^-0.08 


25511° 




80/81 


576i;,J^ 


9 1 q+0.28 
^•l'J-0.24 


57 


-0. 


aq+0.07 

■°^-0.07 


3061^^ 


1 14+0-11 


63/81 


6931?? 


44+0.34 
'^•**-0.30 


58 


-0. 


07+0. 07 
■°'-0.07 


287lfg 


i^-i'^-o.ii 


72/81 


6501^^ 


90+0.31 
^•^•J_0.27 


59 


-0. 


qt5+o.07 

■^^-0.06 


3541^3^ 


-1 97+0.11 
^•^ ' -0.10 


69/81 


80ll?^ 


t-lU-0.33 


60 


-0. 


7r,+0.06 
■ '"-0.06 


517ir, 


l^-^"-0.09 


81/81 


11711^? 


8.69l°:^^ 


61 


-0. 


74+0.06 
■ '^-0.06 


465li 


1 09+0.11 


77/81 


10521?° 


r 79+0.49 
■J- '^-0.45 


62 


-0. 


70+0.07 
■ '"-0.07 


4091^^ 


1 41+0-13 
1^-^1-0.12 


106/81 


9251^^ 


*-l"-0.35 


63 


-0. 


f;f;+0.07 
■'^'^-0.07 


352I26 


1 09+0.12 


56/81 


796l«g« 


4 94+0.38 
^•^^-0.34 


64 


-0. 


09+0. 08 
■0^-0.07 


ooq+37 


1 97+0.13 
1-^ ' -0.12 


63/81 


767l?3 


07+0.32 
'^•"' -0.28 


65 


-0. 


or+0.09 
■°'^-0.08 


270l^^ 


1 19+°" 


72/81 


613l?2 


2 S4+0-27 

^"^^-0.23 


66 


-0. 


71 +0.07 
■ ' -'--0.07 


3561^1 


1 20+°-" 
i-^"-o.io 


73/81 


806l^°3 


4 Q-l+0-39 
*"Jl-0.35 


67 


-0. 


6Q+0-07 
■"^-0.07 


000+30 

000 


1 97+0.12 
l-^'-O.ll 


91/81 


75416? 


Qri+0.32 
O.OU_Q 28 


68 


-0. 


aq+0.07 

■°^-0.07 


3881^^ 


-1 9-1+0.12 
1-^1-0.11 


91/81 

18 , 
65/81 


8771^?° 


9r+0.35 
'3-^J_0.3l 


69 


-1. 


cn+0.30 
■'J"-0.20 




1 ^6+°-3'^ 

1.0U_g 5g 


or -1+294 
^^1-89 


0.38l°:J^ 



f: constant factor of the WAM against the BAT in joint spectral fittings. 
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Fig. 8. The relation between isotropic equivalent luminosity Liso and the source 
frame -E-pcak measured from each 1-scc timc-resolvcd spectra as shown in ta- 
ble 2 and 3. The power law bet fit model is also shown by red dotted line. 
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Fig. 9. The light curve (upper panel) and -E'pj,g^j./L?Q as the function of time (lower panel). Vertical ar- 
rows in top panel represent the outliers from the best-fit power law model in E^^ak'-^lso plot. These 
outliers are located above the 3 a limit of the best-fit constant parameter in the lower panel. As 
seen, all outliers belong to initial rising phase of each pulse. The best-fit constant parameter and 
3 cr confidence level arc also shown by horizontal solid and dashed lines in lower panel, respectively. 
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Fig. 10. Same as Figure 8 but divided into two populations; main popula- 
tion (black circles) and outliers (red squares). The power law best- fit model 
for each population are also shown by black solid and red dashed lines. 



data and folded model 




channel energy (keV) l-.so(''0'^^ 6rg s-') 

Fig. 11. Left panel shows a comparison of the integrated spectrum of all outliers (red squares) with that 
of several main population (black circles). The selected data are indicated by the thick points in the 
Etp — Ljso plane shown in the same figure. The differences of the spectral shape for these two population 
are clearly seen. Right panel shows the relation between the L\^^ and the E'tpeak value obtained by the 
individual simulated spectra (black circles), and by the total integrated spectrum (red square). This figure 
indicates that the measured S'peak value could be higher only by factor about 1.3 even if the hard-to-soft 
evolution following the Yonetoku relation exists at a finer time scale in initial rising phase of each pulse. 
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flux normalized by the time-averaged value Energy (keV) 



Fig. 12. Left panel shows a relation between burst intensity and the spectral parameters ob- 
tained by the Monte Carlo simulation. The burst intensity and obtained parameters are normal- 
ized by the value obtained by the time-averaged spectral analysis. The 20% unceratinty region 
is indicated by the dashed-lines. This figure indicates that the uncertainty of the spectral pa- 
rameter depending on the burst intensity should be within 20%. Right panels shows the com- 
parison of the background spectra between before (To(WAM)-52 s to To(WAM)-2 s) and after 
(To(WAM)-|-87 s to To(WAM)-|-137 s) time interval. The residuals of these spectra are shown in 
the lower window. The change of the background level is less than 5% during this observation. 
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